The Origin of Stokes Shift in Semiconductor Quantum Dots 
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Abstract 

The mechanism of Stokes shift in semiconductor quantum dots is investigated by calculating the energy of the excitonic 
states. We have taken into account all possible contributions to the total electronic energy in the dot, i.e, dielectric mismatch 
between dot and surrounding medium, the effects of finite barrier height and electron-hole exchange interaction. The Stokes 
shift is calculated as a function of radius of dot and compared with experimental data on two different semiconductor based 
quantum dots. These results provide evidence for exchange splitting of excitonic states, as the mechanism of Stokes shift in 
quantum dot 

PACS numbers: 73.22.-f, 73.21. La, 72.80.Ey 
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The red shift of the emission spectra with respect to 
absorption spectra is known as Stokes shift, which is com- 
monly observed in semiconductor quantum dots(QDs) 
and the most important quantity that determines the op- 
tical properties of QDs. However the underlying mecha- 
nism of Stokes shift is still controversial. In case of semi- 
conductor QDs, this shift decreases with the increase in 
radius of the dot and disappears beyond a certain radius. 
The red shift, principally, occurs if either the top of the 
valence band is an optically passive state (say P state) or 
if the electron and the hole are in a triplet state. Absorp- 
tion of a photon from the top of the valence band in such 
cases is not allowed and is possible only from an optically 
active state lying deeper in the valence band. The exci- 
ton, once formed after absorption, can not decay to the 
top of the valence band by direct dipole transition and 
hence is called "dark exciton'—. Deexcitation eventually 
takes place with the help of phonons thus giving rise to 
red shifted photons. 

Two different mechanising^^ have been proposed to 
explain the origin of Stokes shift in semiconductor QDs, 
as shown in Fig. 1. Absorption takes place from the 
S-state which lies deeper than the P-state in the valence 
band to form an exciton in the singlet state of the electron 
and hole. According to the first mechanism, (1) of Fig. 
1, deexcitation takes place into the P state with the help 
of phonons giving rise to a Stokes shift by an amount 
AE$ p representing the difference in energy between the 
exciton states formed with S and P hole states(Ref.2). 
According to the second mechanism the exciton in the 
singlet state first thermalizes into a triplet state from 
where it deexcites with the help of phonons either (a) to 
the deeper lying original S state, (2) of Fig. 1, giving 
rise to a Stokes shift by an amount AE$t representing 
the singlet triplet splitting(Ref. 3), or (b) to the P state 
at the top of the valence band, (3) of Fig. 1, with the 
Stokes shift given by AEsp + AE$t- In addition to this, 
the red shift is also observed if the final excited state has 
a different atomic configuration than the initial ground 
state. Such a red shift, known as Frank Condon shift, 
is commonly observed in molecules and point defects in 
solids and may be relevant in ultra small QDs 5 . 

In this Letter we investigate the origin of Stokes shift in 



III-V(InAs) and Il-VI(CdSe) semiconductor based QDs 
by calculating the energies of excitonic states. We have 
taken into account the dielectric mismatch effects, the ef- 
fect of finite barrier height and the effect of the electron 
hole exchange interaction and then compared the results 
with experimental data. InAs and CdSe are chosen be- 
cause of their (i) different crystal structures, (ii) different 
crystal field effects which is zero in zinc blende struc- 
tures^. e InAs) but significantly large in wurtzite struc- 
ture(CdSe), (iii) different effective masses which is small 
in InAs, but large for CdSe and (iv) different nature of 
bonds. It has been observed experimentalljiii^^ that the 
red shift decreases with radius and disappears beyond 
a certain radius. Essentially, the strong dependence of 
Stokes shift on the radius of the QD, is used to inves- 
tigate the origin of Stokes shift in QDs based on these 
semiconductors(InAs, CdSe). 

Optical transitions from the electron state ip e (r) = 

jo(^r)Y o {0,(f>)\Sa), where \Sa) are the Bloch func- 
tions of the conduction band, are possible only to hole 
states in the valence band which have the S-state as a 
component because J Y Lm YoodU = 6 L ^ 5 m , . The optical 
transition probability from the bottom of the conduction 
band to a hole state in the valence band is given bj£ 

P a p = J drr 2 j Q (^-r) (cR (r)) \{u^(3 |ep| Sa) | , where 

c is the component of the S-state R (r) in the hole wave- 
function, (3 is the spin projection of the removed electron 
in the hole state and are the Bloch wavefunctions of 
the hole states at the valence band top. P a p is zero if 
a / /?. Optical transitions from the bottom of the con- 
duction band to the top of the valence band is not possi- 
ble if either (I) the hole state does not have a component 
of the S-state (i.e if c=0), or (2) the spin projections of 
the removed electron in the hole state and the electron 
in the conduction band are not equal, a ^ f3. Since the 
hole spin is obtained by flipping the spin of the removed 
electron, it follows that optical transitions are not possi- 
ble if the spins of the electron and the hole are parallel 
i.e if they are in a triplet state. Thus, if the top of the 
valence band is a P-state, or the electron and the hole 
are in a triplet state, a dark exciton is formed. 
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The exciton states are obtained from the Hamiltonian 

H = H e + Hh + H so + V e -h + Vpol-s + Vpol-eh (1) 

where H e , H^, H so ,V e -h are the electronic Hamiltonian, 
the hole Hamiltonian, the spin-orbit interaction and the 
e-h coulomb interaction respectively given in Ref.2 and 
7. The last two terms represent the surface polarization 
energies arising due to the difference in the dielectric con- 
stants between the semiconductor quantum dot and the 
surrounding medium-. Vp i- S is the self energy of the 
electron and hole due to their image charges and Vp i-eh 
is the mutual interaction energy between the electron 
and hole via image charges. The excitonic states ob- 
tained from Eq.QJ will be split further due to exchange 
interaction- H ex between the electron and hole. 

Since the nature of the state at the top of the valence 
band plays a pivotal role in the formation of dark and 
bright excitons, we first investigate the hole states for 
both the zinc-blende and wurtzite structures using the 
k.p hole Hamiltonian Hh + H so given in Ref.9. The en- 
ergy eigenvalues of the hole states without inclusion of 
excitonic effects as a function of dot radius for InAs (zinc- 
blende structures) and CdSe(wurtzite structures) are 
shown in Figs. 2(a) and 3(a) respectively. These eigen- 
values get modified when coulomb interaction between 
electron in the conduction band and hole in the valence 
band are taken into account. The corresponding results 
including excitonic effects are shown in Figs. 2(b) and 
3(b) respectively. 

From Fig. 2(a) we observe that in InAs the optically 
passive Pa is the ground state which lies below the first 
optically active state S| state. The states are labeled 
Sa, P3 etc where the capital letters correspond to the 
lowest L present and the subscripts gives the total an- 
gular momentum F'=L + I + S, where L is the orbital 
angular momentum, I(=l) is angular momentum of the 
Bloch wavefunction at the valence band top and S is the 
spin of the hole. The difference in energy between the op- 
tically passive P| state and the optically active Sa state 
decreases as the radius of the dot increases and eventu- 
ally S3 becomes the ground state at around R ~ 50A. 
Hence dark exciton formation and associated Stokes shift 
will be observed below this radius. 

For the wurtzite structure, the Hamiltonian is axially 
symmetric and only the z-component M = L z + I z + S z of 
the total angular momentum is conserved. From Fig. 3(a) 
it is observed that for CdSe the optically passive \P X t) 
with M = \ state is the ground state at R ~ IOA. The 

energy of \P X f) in units of e c = 2 (71 is a Lut- 

tinger parameter) has an upward slope as a function of 
radius and crosses the relatively flat \S X T) state with 
M = I at R r-j 28A. Dark exciton formation and Stokes 
shift is predicted for CdSe below R ~ 28A. The trend of 
S-state and P state can be explained by observing that 
\S X T) with M = I is predominantly made from only one 
state with composition L = 0, I z = +1,S S = +§ whose en- 
ergy decreases as but in units of e is independent 



of R and hence the state is relatively flat. But the op- 
tically passive P state with M = | can be made either 
from the state with composition L z = —l,I z = +l,S z = 
+I(j 2 = i z + s z = +§) or the state with the composition 
L z = +i,h = —1,S Z = +\{Jz = -5). These two states are 
separated from each other by spin-orbit interaction. At 
very low radii spin orbit effects are negligible relative to 
kinetic energy terms due to confinement. Hence the two 
states are almost degenerate at R ~ 10A. The coupling 
between these two states give rise to a splitting of the 
states with the lower P state lying below the S state. 
As the radius increases the energies of both the S and P 
states decrease. The energy of S state decreases as 
but the energy of the P state decreases more slowly be- 
cause of the presence of the spin-orbit interaction energy 
which becomes important as the kinetic energy starts de- 
creasing with the increase of radius. Eventually the two 
levels cross. In Figs. 2 and 3 the energies have been 
shown in the units of e Q . 

When the coulomb interaction between the electron in 
lS e state and the hole is taken into account, the energies 
decrease. The decrease in energy of S-state is more than 
that of P state, the overlap between lS e state and S 
state of the hole being greater. This causes the crossing 
between S and P states to occur at lower values of R as 
can be seen from comparison of Figs. 2(b) and 3(b) with 
the corresponding Figs. 2(a) and 3(a). The differences 
between the excitonic S and P states AEgp as a function 
of R, are shown in the insets to Figs. 4 and 5. 

The surface polarization charges arising due to the di- 
electric mismatch between the dot and the surroundings 
also affect the excitonic energies. The contributions of 
self energy Vp i- S and attractive energy Vp i- e h are op- 
posite in sign and comparable in magnitude^. Hence the 
dielectric mismatch effects are greatly reduced with the 
net effect being an increase in energy depending on the 
amount of dielectric mismatch between the dot and the 
surrounding. The increase in energy for the P state is 
more than that for the S state resulting in a crossover 
to occur at a lower dot radius. This can be seen from 
Figs. 2(b) and 3(b) which gives the comparison between 
zero and finite dielectric mismatch, i.e e = 1 and e = 3. 
The self energy of a charge confined to a spherical dielec- 
tric is minimum—i at r=0. Since (r 2 ^ for the P-state is 
larger compared to the S-state, the self energy Vp i- S is 
higher for the P-state. But the attractive energy Vp i- e h 
is higher for the S state compared to the P state. This 
leads to a higher net increase of polarization contribution 
[Vpoi-s + Vp i- e h] for the P-state as compared to the S- 
state and hence to a shifting of the S-P crossover to lower 
values of R as can be seen from Figs. 2(b)and 3(b). 

Another cause of dark exciton formation and the con- 
sequent Stokes shift is the splitting of the active ex- 
citon states due to electron hole exchange interaction 3 

Hex = (-§) e ex (a Q ) 3 8 (r e - ffA cr.J where t ex is the ex- 
change strength constant, a the lattice constant, a the 
Pauli spin matrices representing the electron spin and J 
is the hole spin matrix. The exchange interaction gives 
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rise to a splitting between triplet and singlet state of the 
electron and hole, with the triplet state lying lower in 
energy. The triplet state being passive, excitation takes 
place to the higher singlet state. Deexcitation can take 
place from the triplet state reached by thermalization, 
thus giving rise to a red shift of the emmision spectra 
with respect to the absorption spectra. 

In InAs the first optically active S state (Fig. 2(a)) has 
total anguar momentum F' = |(L = 0, 1=1, S=^) and 
is fourfold degenerate. The lS e electron state in the con- 
duction band being doubly degenerate, the exciton states 
are eightfold degenerate in the absence of exchange inter- 
action. The exchange interaction splits these degenerate 
levels into two: the lower optically passive triplet state 
with F=2(F = F' + S e ) and the upper optically active 
singlet state with F=l (Fig. 2(c)). The difference in 
energies AEst between the singlet and triplet states as 
a function of dot radius is shown in Fig. 4. AEst de- 
creases with R and vanishes at R ~ 35A. The exchange 
interaction is enhanced with decrease in radius because 
of the increase in the overlap of electron and hole wave- 
functions due to confinement. 

In CdSe (Fig. 3(c)) the ground state with F z = ±2 can 
be obtained by coupling a hole with M = ±| with the 
electron state with (S e ) z = ±|. Since these combinations 
form triplet states, the ground state with F z = ±2 is 
optically passive and forms a dark exciton. The first 
excited state with F Z =±1 L can be obtained either by 
coupling M=±|, (<S e )z = T-| or alternatively by M=±| 
with (5' e ) z =±l. The state ±1 L will be a mixture of 
singlet and triplet states and will be optically active. The 
difference in energies AEst between the optically passive 
ground state and the first optically active excited state 
as function of R for CdSe is shown in Fig. 5. As in the 
case of InAs, AEst for CdSe also decreases with R and 
vanishes around R ~ 50A. 

As seen from Figs. 2, 3, 4 and 5 the curves correspond- 
ing to finite barrier case lie lower than infinite barrier 
case. Both the exchange interaction and coulomb inter- 
action depend on the amount of overlap of electron and 



hole wavefunctions, which is reduced in the finite bar- 
rier case due to penetration of the substantial amount 
of electron wavefunction outside the dot resulting in the 
decrease of exchange energy and coulomb binding energy. 

Having theoretically calculated AEsp and AEst we 
are in a position to investigate the various mechanisms 
for Stokes shift given in Fig. 1 by comparison with ex- 
perimental data. In Figs. 4 and 5 are shown the plots of 
AEst as a function of the dot radius together with the 
experimentally observed Stokes shift for CdSe and InAs 
respectively. AE$p as a function of R for these cases 
are shown in the insets of the corresponding figures. It is 
observed that AEs p is much larger than the experimen- 
tal red shift for InAs. For CdSe, as seen from the inset 
of Fig. 5 at R ~ 10A AEsp is almost the same as the 
experimental Stokes shift but for larger values of R the 
experimental values are much higher. It is seen from Figs. 
4 and 5 that AEst fits the experimental data rather well. 
We observed that the fit for CdSe corresponds to an in- 
finite barrier at the dot radius but for InAs a good fit 
is obtained for a finite barrier height. This can be ex- 
plained by observing that the effective mass of electron 
in InAs is very small and hence the lowest eigenvalue lies 
near the top of the barrier height and there is substantial 
probability of finding the electron outside the dot. For 
CdSe, the effective mass of electron is large and the low- 
est eigenvalue lies near the bottom of the potential well, 
making the infinite barrier effects predominant. In Fig. 5 
we have also plotted the theoretical curve for Stokes shift 
obtained by Efros et a\m in which the coulomb interac- 
tion between the electron and hole was not taken into 
account. We observe that better agreement with exper- 
imental Stokes shift is obtained by our theoretical curve 
which takes into account the coulomb interaction. 

From a comparision of the theoretical results with ex- 
perimental data we conclude that the Stokes shift ob- 
served in QDs is caused by the splitting of the exciton 
states by the electron hole exchange interaction through 
the deexcitation mechanism (2) in Fig. 1. 
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Figure Captions 

Fig.l: The schematic representation of excitation- 
deexcitation processes taking place at the band edges. 
Fig. 2: The energy levels in units of £o= {t>) f° r zmc 
blende InAs QDs as a function of R in the (a) absence 
and (b) presence of e-h coulomb interaction(without ex- 
change). The energy reference is taken as the top of 
the valence band without the spin orbit effects included. 
Eigenvalues in the absence(e = 1) as well as presence(e = 
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3 of dielectric mismatch effects are also shown in (b) 
for both the cases of infinite and finite barrier where 
e _ — tdot — _ ^ eigenvalues including c-h coulomb 

and exchange interactions. The states are labelled with 
F z with superscripts standing for upper and lower state. 
The midpoint of the exciton energies with M=±| and 
M=±5 is taken as the reference energy (zero of energy). 
Fig. 3: The energy levels in units of e G for wurtzite CdSe 
QDs as a function of dot radius R in the (a) absence 
and (b) presence of e-h coulomb interaction and dielectric 
mismatch effects (without exchange) for both the cases of 
infinite and finite barrier. The energy reference is taken 
as the top of the valence band with the spin orbit effects 
included. The states are labelled as S x 1, P x t etc, where 
capital letters represent the dominant L present and the 



subscripts stand for X-, Y-like or Z-like states, the arrow 
indicating the spin state, (c) eigenvalues including the 
e-h coulomb and exchange interactions. 
Fig. 4: Singlet triplet splitting in InAs in the absence 
e = l and presence e — 3 of dielectric mismatch effects 
for both the cases of infinite and finite barrier. The inset 
gives the energy differences of S and P states. Exper- 
imental data(empty circle and square) are taken from 
Rcf.4. 

Fig. 5: The singlet triplet splitting in CdSe in the absence 
and presence of dielectric mismatch effects for both the 
cases of infinite and finite barrier. The inset gives the 
energy difference of S and P states with M = | and | 
respectively. Experimental data(empty circle) are taken 
from Rcf.3. 



4 



Exciton 
states 



Valence 

band 

states 



(3) 



(2) 



S singlet 
S triplet 



(1) 



P state 



S state 



Figure 1 




-10 1 

20 „ 40 

Radius (A) 



Figure 2 




-60 1 1 

20 , 40 60 

Radius (A) 



Figure 3 




o i 

10 15 o 20 25 

Radius (A) 



Figure 4 




Radius (A) 



Figure 5 



